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We report on searches for Flavor-Changing Neutral Current (FCNC) decays in pp collisions 
at a center-of-mass energy of 1.96 TeV using up to 0.78 fb _1 of data, collected at the CDF 
' and D0 detectors respectively. The rare FCNC decays presented here are the searches for 

r S^ t B'ld ~* and B° — * (f> fi + n~ . 




1 Introduction 



In the Standard Model (SM), Flavor-Changing Neutral Current (FCNC) decays are highly 
suppressed and can only occur through higher order diagrams. For example the branching 
fractions of the FCNC decays B° sd — ► y + y~~ have a SM expectation^ of B{B® — > y + y~) = 
(3.42 ± 0.54) x 1(T 9 and B{B° d -^'y+y~) = (1.00 ± 0.14) x 10~ 10 . These predictions are about 
two orders of magnitude smaller than the current experimental sensitivity. However, new physics 
contributions can significantly enhance these branching fractions. In the absence of an observa- 
tion, any improvements to the limits can be used to set significant constraints on various models 
beyond the SM. The exclusive FCNC decay B® — > 4> y + y~ is related to the quark-level transition 
of b — ► s £ + £~ . The long term goal is study kinematic properties for these decay like the invariant 
di-muon mass distribution or the forward backward asymmetry. But also an observation of this 
decay or experimental upper limit on its rate will yield additional important information on the 
flavor dynamics of FCNC decays. 

2 Search for the rare decays B® d — > y + y~ 

The most stringent published limits are a result of a combination of the recent D0 results using 
0.3 fb _1 of data and the latest published limit from cdfEI using 0.36 fb 1 of recorded data. The 



obtained limits at 95% C.L. were B{B® -» < 1.5 x 10~ 7 and -» < 4.0 x 10~ 8 

respectively. 

The used analysis procedure is similar in the D0 and CDF experiments. After a pre-selection 
step, special discriminating variables in an optimization procedure are used to further enhance 
the signal efficiency while reducing the expected background. Both collaborations employ a 
"blind" analysis strategy when choosing the final selection criteria. As normalization the well 
measured and high statistics decay B + — > J /ipK + was used. In the following the most recent 
analyses from both experiments will be described. 

The strategy of the D0 collaboration is to use the already published analysis procedure^ for 
the data set recorded after the first publication, to obtain a sensitivity. This obtained sensitivity 
is then combined with the already obtained limit. CDF used a published analysis procedure^ 
as well and added more data to get an updated limit. 

D0 used data collected by dimuon tiggers with two muons of opposite charge that form a 
common secondary 3D-vertex with an invariant mass between 4.5 and 7.0 GeV/c 2 . Each muon 
candidate had to have px > 2.5 GeV/c, \t]\ < 2.0 and a sufficient number of hits in the central 
tracking station. To ensure a similar pp dependence of the p + p.~ system in the signal and in 
the normalization channel, p^ had to be greater than 5 GeV/c. 

For the D0 analysis the following variables were found to exploit best the distinctive proper- 
ties of the B meson: the 3D Opening- Angle, A0 (here the decay length vector is determined in 
three-dimensions), the 5-candidate Isolation, /, and the Transverse Decay Length Significance, 
L2D / &l 2 d- After optimizing the cuts and a linear interpolation of the sideband population for 
the data sample into the signal region we obtain an expected number of background events of 
2.2±0.7. The remaining background distribution is shown in Fig. ^ (left side). The values of 
the discriminating variables for this additional data set changed only slightly with respect to 
the previous analysis. The resulting sensitivity is calculated using the B + — > J/i/jK + decay as 
normalization. The invariant mass distribution of the decay channel is shown in Fig. ^ (right 
side). This yields a sensitivity including all the statistical and systematical uncertainties of 
(B(B® -> n+fi')) < 3.0(3.7) x 10- 7 at the 90% (95%) C.L. 
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Figure 1: The invariant mass distribution for the decay B® — > fi + fj,~ after optimized requirements on the discrim- 
inating variables on the left side. On the right side the invariant mass distribution of the normalization channel 

B + — » J/ipK + is shown. 

The obtained expected limit for the new data set is combined with the previous analysis^ 
using 0.3 fb^ 1 of data. The expected limit for the entire data set is thus a combination of two 
measurements: one actual measurement with observed events (old data) and one hypothetical 
experiment (new data) with all possible outcomes in the signal region weighted by their Poisson 
probability of occurrence, i.e., under the assumption of a background hypothesis only. We 
have used a Bayesian technique to combine the two experiments and include at this stage the 
uncorrelated uncertainties into the expected limit calculation. The background uncertainty, the 
uncertainty on the number of normalization events and the statistical error on the efficiencies 
are considered as uncorrelated. The expected upper limit at the 90% (95%) C.L. for the entire 



D0 data set of 0.7 fb~ 
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-> < 1.9 (2.3) x 10 
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The CDF detector has of two muon detectors with different acceptances. Four layers of 
planar drift chambers (CMU) detect muon candidates with px > 1.4 GeV/c and provide coverage 
in the pseudorapidity range \rj\ < 0.6. The central muon extension (CMX) consists of conical 
sections of drift tubes and extends the coverage to 0.6 < \rj\ < 1.0 for muon candidates with 
Pt > 2.0 GeV/c. The data used in the analysis are selected by two classes of dimuon triggers: for 
the CMU-CMU (U-U) triggers both muon candidates are reconstructed in the CMU chambers, 
while for the CMU-CMX (U-X) triggers one of the muon candidates is reconstructed in the 
CMX chambers. Since they have different sensitivities, the U-U and U-X channels are treated 
separately, combining the results at the end. The offline reconstruction begins by identifying two 
muon candidates of opposite charge which satisfy the online dimuon trigger requirements and 
with pt > 2 (2.2) GeV/c. The random combinatoric backgrounds are suppressed by requiring 
the vector sum of the muon transverse momenta to be \p^\ > 4 GeV/c. The 3D decay length 
is given by L^d = L ■ p^ /\p^\, where L is the displacement vector from the primary to the 
dimuon vertex. The primary vertex is determined using a constrained vertex fit of all tracks in 
the event, excluding the n + pair and other secondary decay tracks. For each £>-candidate the 
proper decay time r = M^^L^/lp^l is estimated, where M w is the invariant mass and p^ 1 is 
the momentum vector of the dimuon system. Additional background is reduced by demanding 
L3B < 10 cm, the uncertainty on L3D to be less than 150 //m, and 2o~\ < A < 0.3 cm, where 
A = c • r and o~\ is the total uncertainty on A. 

To enhance signal and background separation CDF constructs a multivariate likelihood ratio 
based on the input variables: I, AO, and A probability P(A) = e X ^ CTB "W ^ w here Tb is the 
world average B s ^ lifetime. CDF uses the P(X) variable instead of A in constructing the 
likelihood ratio because the -P(A) distribution is nearly flat, and therefore better behaved in the 
likelihood. The likelihood ratio is then defined to be 



where x\ = I, x<z = AO, x% = -P(A), and P s ^(xi) is the probability that a signal (background) 
event has an observed X{. The probability distributions for the signal events are obtained from 
the signal MC and the background distributions are taken from the data sidebands. An example 
of the likelihood ratio for the U-U channel is shown in Fig [2] (left side) . The optimization of 
the analysis is based on the a priori expected 90% C.L. upper limit on B(B® d — > fi + /i~). The 
expected limit for a given set of optimization requirements is computed by summing the 90% C.L. 
limits over all possible observations N Q , weighted by the corresponding Poisson probability when 
expecting A5 background events. A scan over a range of Lr requirements the optimal value was 
determined to be Lr > 0.99. With the optimized selection requirements, the expected number 
of background events in 0.78 fb _1 for B® in a signal region of ±60 MeV/c 2 (±2.5o"m) around the 
world average B° mass 5 is N b is 0.88±0.30 [0.39±0.21] for the U-U [U-X] channel. The expected 
number of background events in a ±60 MeV/c 2 signal region around the B® is estimated to be 
N b is 1.86 ± 0.34 [0.59 ± 0.21] in the U-U [U-X] channel. Using these criteria 1 [0] and 2 [0] 
events are observed in the U-U [U-X] channel of the B® and B® signal boxes respectively which 
is consistent with the background expectations, shown in Fig. [21 (right side). Combining the 
U-U and U-X channels taking into account the correlated uncertainties, 90% (95%) C.L. limits 
can be derived to be 
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B(B° S -> n + n~) < 8 x 10~ 8 (1.0 x 10~ 7 ) and 
B{B° d -> /iV~) < 2.3 x 10~ 8 (3.0 x 10~ 8 ). 
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Figure 2: Signal and background distribution of the likelihood ratio for the U-U channel on the left side. Invariant 
mass distribution versus the event likelihood with the signal boxes indicated on the right side. 



3 Search for the rare decays B® — > 4>fJ, + fJ- 



The decay Bg — > /U + /x is an exclusive FCNC in the B® meson system. Within the SM, 
the decay rate for the B® — > (ft fi + fi~ decay, neglecting the interference effects with the much 
stronger 5° -> JM and B° s -» ifi(2S)(f) resonance decays, is predicted to be of the order of 
1.6 x 10~ 6 (see Ref.^f) with about 30% uncertainty due to poorly known form factors. 

A search for this rare decay is performed by the D0 collaboration and uses about 0.45 fb _1 of 
data. The event pre-selection starts by requiring exactly two muons fulfilling quality cuts on the 
number of hits in the muon system and the two additional charged particle tracks to form a good 
vertex. The invariant mass of the two muons is required to be within 0.5 < vn^^- < 4.4 GeV/c 2 . 
In this mass region, the J/ip(-^> fi + 'fi~) and iJj(2S)(—> resonances are excluded to dis- 

criminate against dominant resonant decays. The x 2 /d-o.f. of the two-muon vertex is required 
to be less than 10. The tracks that are matched to each muon are required to have a sufficient 
number of hits in the tracking detectors and the transverse momentum of each of the muons 
is required to be greater than 2.5 GeV/c with \rj\ < 2.0 to be well inside the fiducial tracking 
and muon detector acceptances. The number of B® — > <f) A t+ / U_ candidates is further reduced 
by requiring pj, > 5 GeV/c and asking the B® candidate vertex to have x 2 < 36 with 5 d.o.f. 
The two tracks forming the cp candidate are further required to have px > 0.7 GeV/c and their 
invariant mass within the range 1.008 < < 1.032 GeV/c 2 . In analogy to the search for 
the decay B® — > // + /i~the same discriminating variables and optimization procedure is used to 
enhance the sensitivity. The optimal values of the discriminating variables yield after a linear 
interpolation of the sidebands into the mass window signal region, 1.6±0.4 expected background 
events as shown in Fig. [3] (left side). Upon examining the data in the mass region, zero candi- 
date events are observed in the signal region, consistent with the background events as estimated 
from sidebands. The Poisson probability of observing zero events for an expected background 
of 1.6 ± 0.4 is p = 0.22. In the absence of an apparent signal, a limit on the branching fraction 
B(B® —s- (j)fi + fi~) can be computed by normalizing the upper limit on the number of events in 
the B® signal region to the number of reconstructed B® — ► J/ipcf) events. The mass spectrum of 
the reconstructed B® — > J/ip<i> after all cuts is shown in Fig. (right side). Including statisti- 
cal and systematic uncertainties, without including the uncertainty on the measured branching 



fraction of the decay B® 



■ J/if)(f>, the Feldman and Cousins (FC) limit is 



< 4.4(3.5) x 10" 



at the 95% (90%) C.L. respectively. Taking a Bayesian approach with a flat prior and the 
uncertainties treated as Gaussian distributions in the integration, an upper limit of B(B® — > 
0H + fi-)/13(B° -> J/i/}</>) < 7.4(5.6) x 10" 3 at the 95% (90%) C.L., respectively was found. 
Using only the central value of the world average branching fraction^ of B(B® — > J/ifxfi) = 
(9.3 ± 3.3) x 10" 4 , the FC limit corresponds to B(B° S -> <f> fi + /J,') < 4.1 (3.2) x 10~ 6 at the 95% 
(90%) C.L. respectively. This is presently the most stringent upper bound and has been recently 
submitted to PRlM. 
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Figure 3: The invariant mass distribution for the decay B® — > /i" 1 "/^ on the left side. On the right side the 
invariant mass distribution of the normalization channel B® — > J/ip4> ls shown. 
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